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PACS 98.62.Mw - Infall, accretion, and accretion disks 

PACS 52.35.Py - Macroinstabilities(hydromagnetic, e.g., kink, fire-hose, mirror, ballooning, tear- 
ing, trapped-particle, flute, Ray leigh- Taylor, etc. 
PACS 47 . 65 . Md - Plasma dynamos 

Abstract. - Differentially rotating plasma disks, around compact objects, that are imbedded in a 
"seed" magnetic field are shown to develop vertically localized ballooning modes that are driven 
by the combined radial gradient of the rotation frequency and vertical gradients of the plasma 
density and temperature. When the electron mean free path is shorter than the disk height and 
the relevant thermal conductivity can be neglected, the vertical particle flows produced by of 
these modes have the effect to drive the density and temperature profiles toward the "adiabatic 
condition" where t)t = (dlnT/dz) / (dlnn/dz) — 2/3. Here T is the plasma temperature and n the 
particle density. The faster growth rates correspond to steeper temperature profiles (vt > 2/3) 
such as those produced by an internal (e.g., viscous) heating process. In the end, ballooning 
modes excited for various values of t)t can lead to the evolution of the disk into a different current 
carrying configuration such as a sequence of plasma rings. 
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I. Introduction. — Axisymmetric plasma disks [1,2] 
in the prevalent gravity of a central object that do not have 
internal currents and are imbedded in a "seed" vertical 
magnetic field can be described in their equilibrium state 
by one-fluid equations that are basically the same as those 
describing gaseous disks [3]. Then axisymmetric plasma 
ballooning modes that are contained vertically within the 
disk and are radially modulated over the "Alfven scale 
distance" (v^/fi) are identified as the most likely process 
leading the starting disk configuration to evolve into a 
current carrying one such as the axisymmetric sequence 
of rings found in Ref . [4] . The driving factors of the fastest 
modes, within this category of modes, are the combined 
effects of the vertical particle density and temperature gra- 
dients and the radial gradient of the rotation frequency. 
The particle motion is assumed to be "frozen" to the mag- 
netic field. 



The relevant quasi linear theory shows that this kind 
of mode, that is referred to as the "thcrmo-rotational 
instability" , produces a vertical particle inflow toward 
the equatorial plane and a temperature outflow if = 
(dlnT/dz) /(dlnn/dz) > 2/3. The reverse occurs if t\t < 
2/3 . Here T is the plasma temperature (2T = p/n), p 



is the total plasma pressure and n is the particle density. 
Therefore, the main effect of these modes is to drive the 
density and temperature profiles toward the "adiabatic" 
condition t\t = 2/3. When this condition is reached the 
slowly growing ballooning modes presented in Ref. [5] can 
take over and produce, at saturation, a new stationary 
axisymmetric configuration [4]. 

In this context we note that the well known Magneto- 
Rotational Instability (M.R.I.) , driven simply by the gra- 
dient of the rotation frequency Q, was derived originally 
for a cylindrical geometry [6-8] where it can develop a 
growth rate close to fl and produce a significant rate of 
radial transport of angular momentum. On the other hand 
this instability is severely depressed when a strong bend- 
ing of the magnetic field lines, imposed by the geometry of 
a thin disk, is considered. When rjx > 2/3 , correspond- 
ing to a temperature profile for a disk subject to a strong 
heating source (e.g. viscous) around the equatorial plane, 
the growth rate of the ballooning thermo-rotational insta- 
bility is not too inferior to that of the "cylindrical" M.R.I, 
for equal value of the rotation frequency. We note also 
that the thcrmo-rotational instability would not be found 
if the flow compressibility were neglected. 



II. Initial (Stationary) Configuration. We con- 
sider a current-less plasma disk whose density profile near 
the equatorial plane at the reference radial distance Ro 
from the axis of symmetry is represented by 



where 



p[z 



(1) 



where t]t can be positive as well as nul and negative and 
Hq <C Rq- This approximation is consistent with the 
analysis of ballooning modes of the type introduced in 
Refs. [5, 9] that are localized over a height of the order of 
A z such that A 2 Z < Hq. For the sake of simplicity, we 
assume that the disk is immersed in a vertical magnetic 
field and that no appreciable internal currents are present. 
Then the initial state is simply described by the vertical 
and radial components of the total momentum conserva- 
tion that, for a thin disk in which no poloidal velocities 
are present, reduce to 







and 



n (r) = n k 



(2) 



(3) 



where = fl (R) R is the plasma toroidal velocity, 

1/2 

ilk = (GM*/Rq) " is the Keplcrian frequency, M* 
is the mass of the central object, p = 2nT, p = nrrii 
and TOj is the ion mass. In particular, Eq.(2) yields 
Hi = 2To(l + f lT )/(m t n 2 k ). 



7o<U = V0 + Rjftt,R, 



(10) 



and is the perturbed toroidal velocity. Then we follow 
the same proceedure as that in Ref. [5] considering high—/? 
regimes where v\ < <? s , c s is the sound velocity and = 

1/2 

B / (47rp) ' is the Alfven velocity. As pointed out in Ref. 
[5, 9]the present problem involves prominently the scale 
distance V/t/fii,. 

III. Importance of Thermal Fluctuations. The 

z-component of the total momentum conservation equa- 
tion can be written as 



dp 

dz 



= ~57 - W I 11 ) 

and we note that the two terms on the r.h.s. scale as 

P /(pA z ) and pA z /(pH 2 ), (12) 

relative to each other. Therefore, if we consider cases, as 
was done in Refs. [5], where 

(13) 



-zfli, can be neglected 
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the term \pg z \, where g z 
relatively to\dp/dz\ . 



We consider axisymmctric perturbations from the equi- 
librium state represented by 



vj?, = Vr (R - Ro, z) exp (70*) . 



(4) 



and we make use of the following total momentum conser- 
vation equation 




--7-B ■ VB + zptiie z = 0. (5) 



We adopt the "hyperconductivity" condition 
E+(vxB + vxb) /c = 



and write 



(6) 
(7) 



where v p is the poloidal velocity and £ p the relevant dis- 
placement. Thus taking the curl of Eq.(6) we obtain 



d - - d 

Br = B z —£ r , B < f, = B z —t 

az oz 



B z = B z ^-l z -B z ( 



(8) 
(9) 



On the other hand, if \p/p\ 
that 



{ h q/ A D \p/p\ implying 



p/p~-T/T, (14) 
the pg z term has to be retained. Since p/p = p/p + T/T, 

d 
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(15) 



Thus the g z term remains associated with the tempera- 
ture perturbations. 

In this connection we note that a simplified thermal 
energy balance equation that includes an isotropic thermal 
conductivity is 



(70 + D th k 2 R ) f + 7o e 



dT 

dz 



-370 



P rt. 2 dpT 
n dz p 



(16) 



where Z? t ^ is the relevant thermal diffusion coefficient, con- 
sidering perturbations that are appropriate for the modes 
we analyze, as shown in Ref. [5] are locally periodic in 
(R — Ro) and represented by 



f~f(z) exp [ 7o i + ik R (R-Ro) 



(17) 



where k R R,Q ^> 1 and f (z) is localized over A z . Then and, consequently, to 



~p - pVL\z ~ ^-p - zQ 2 k p£ z 



370 1 dp 



57o 



(18) 



where v t h = Dthk R - Given Eq. (11), it is clear from the 
present result that the analysis given in Ref. [5] is valid 
for 



70 > n 2 



3A 



H 2 (5 7 o + 3vth) 



(19) 



considering that the growth rate that is found in this case 
is of the order of 70 ~ Va/Ho- Clearly the inequality 
(19) holds either for relatively large thermal conduc- 
tivities (jo/fth <1) or "adiabatic profiles" such that 
\t]t - 2/3| < 1. 



IV. Ballooning Mode Equations. We define 



r - 6 f 2 



(20) 



and consider it to have significant values of either signs. 
Clearly, flat temperature profiles correspond to Ct — 
—4/5. Conversely, a heating process localized around the 
equatorial plane can be envisioned to maintain relatively 
peaked temperature profiles corresponding to Ct > 0. We 
refer to Eq. (5) and note that the displacement £^ is re- 
lated to £r by the <p - component of the total momentum 
conservation equation that is 



1 (7o6 



20*6 



dz 2 ^ 



We consider the limit where 



(21) 



(22) 



that will be justified a posteriori. Therefore Eq.(21) re- 
duces to 

d 2 ~ 

2 l0 n k £ R ~v 2 A — (23) 

The d (e^, ■ V x A m ^ /dz = equation derived from the 

total momentum conservation equation A TO = can be 
written as 



dz 2 
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We reconsider Eq. (14) that leads to 

V ■ £ ~ - (3/5) L (dp/dz) /p ~ I A z /H 2 



= 0. (24) 



(25) 



9 c „ 



(26) 



Then, making use of Eq. (23) we obtain 
p (j*£ R + 2Sl k n' k R £ R 



dz 2 



4tt V dz 2 



k 2 



P7o(4 + ^Cr-^hO^(. 



' dz 



where 



Therefore, the r.h.s. of Eq. (27) becomes 



'H 2 




where 



dz \ k R I 



(28) 



(29) 



(30) 



Now we define A z = i£, z /k R considering the case where 



kW 



R A 



-2n k n' k R a - 2fc \5k R \v 2 A , k > 0, for k 2 



-2Qfcf2j.i?o/v^, and localized modes for which p can be 
represented by Eq. (1), Eq. (27) reduces to 



dz 2 



3Q 2 — 



(2k \Sk R \v 2 A ) 



( ^7o 2 ^ 



(31) 



3f2 2 . By using Eq. (30) this can be 



as — 2R a n' k il k 

rewritten as a quintic equation for A z that can be inte- 
grated to give the quartic 



d 7 



1 d 3 - 



dz\[ ek H 2 ) dz Az+ k 2 dz* A '' 



7 2l 

3v 2 A 



C 7 



'3H 2 



A z = 0, 



(32) 



where = 2 \Sk R \ /ko. The related quadratic form that 
can be used to evaluate 7q is 




(33) 



where () = J H ° Ho dz, we have considered modes that are 
localized over A z < Hq and carried out appropriate inte- 
grations by part taking into account that A*dA z /dz — > 
and A*d 3 A z /dz — > as z 2 — > iJp. Clearly, 7q is always 



real, that is, either purely growing or stable oscillatory 
modes can be found in the context of the adopted approx- 
imations. 

The relevant expression for p (z) is obtained from the 
equation 



(3eo) 2 ^ 3 , where e = 1/ (k H ) < 1. In particular A 2 
(H /k 2 ) 1/3 3 1 / 6 and Eq. (38) can be re-written as 

d 2 



d 2 ~ 

(Jrp Al 

T dk 2 



(3e ) 



2/3 



dk 2 



(kA k ) 



Po 7o + 



dz / ko ' 



derived from Eq. (11) that gives p as a function of £ 2 . 

V. Fourier Transform. — The order of the equation 
to be solved can be reduced further, as we look for verti- 
cally localized modes, by taking the Fourier transform of 
Eq. (32). For this, we adopt the following dimensionless 
variables 



z 

ai 



^(fc A 2 ) 4 , E° k ^e k (k A z ) 2 



(35) 



where (koA z ) 3> 1, and rewrite Eq. (32) as 



d_ 

dz 
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A Z =Q. 



Then we define 



A 



2 _ Hq 



_ r> (^A 2 



k 



and C!p = Ct 



A, 
A 



(36) 



(37) 



and obtain the following equation for the Fourier trans- 
form A k of A z 



d 2 ~ 

r° — A, 

CT dk 2Ak 



A 



dk- 



(kA k ) 



+ [T 2 - k 2 {El - k 2 )] A k 
The relevant quadratic form 



(38) 




(39) 



may be more convienient than Eq. (33) in order to eval- 
uate Tg, for given values of E® and C^., using appropri- 
ate trial functions. Here () = J _ dk, and in deriving 
Eq. (39) we have considered that, for the localized solu- 
tions of Eq. (38), A k dA k /dk and kA k d (jtAkj /dk — > as 
k 2 — > oo. 

Then, if we refer to cases for which \Cr\ ~ 1 wc 
may choose A 2 = 3A 2 /fco- Consequently, (A z /Aq) = 



[(rg - Elk 2 ) + ¥] A k = o. 



(40) 



If instead, we consider cases where |CV| < 1 we may choose 
A 2 = A . Then C£ = C T j (3e ) and Eq. (38) becomes 



r o d 2 A k 



k£-[kA k 
dk 2 1 



[T 2 -k 2 (El-k 2 )]A k = 0. 



(41) 



The cases where IC^I 



1 correspond to the "ordering 
of maximum information" where all the terms in Eq. (41) 
are of the same order. Consequently, these are the cases 
that we shall analyze most closely. We note also that 



€r = *kA k 
and that, referring to Eq. (21), 



To +k 



2 vi 



-i2Vl k ^kA k . 



(42) 



(43) 



Therefore, when we look for even solutions of Eq. (38) 
with A k = A1 ^ for k = it is important to include 
the term in Eq. (43) in order to avoid a singularity 
for Q; at k = 0. This means that for these cases £^ will 
be localized over scale distances of the order of Va/Jo 
and when these become of the order of H , the assumed 
approximation (1) for p is no longer valid. 

VI. Peaked Temperature Profiles. At first we 
analyze the case where the temperature profile is relatively 
peaked and corresponds to r\T > 2/3. We refer to Eq. (41) 
and consider the lowest order solution 



Then we have 

1)- 
which gives 



Crp (k 



r 2 



kik 6 



A° k exp | 



-J? /2) 



(44) 



3k) 



k 2 E° k + fc 4 = (45) 



2 V3v A 2 
70 i^5 (3??T 



and 



El 



a 



Ti 



( k = 3e 



Vt 



2)Q 



(46) 



(47) 



Clearly, this mode is unstable only if t)t > 2/3. We note 
that, the marginally stable (T = 0) modes that can be 



found for = are given by the simple equation 



dk 2 



kA k \- El-k 2 



(kA k ) = 0, 



(48) 



T..A 



which is the Fourier transform of that given originally in 
Ref. [9], 



that gives 



,2 d^ R 
"dz 2 " 



e k U = 0. 



(49) 



Clearly, when Cj. < the considered mode is purely 
oscillatory. When i]t ~ 1 the growth rate that can be 
obtained from Eq. (41) for Cj. ^> 1 and E® 1 can be 
considerable and not much smaller than that of the MRI 
instability in a cylindrical plasma [6]. We note also that, 
for this mode, ^ as given by Eq. (43) requires that the 
term 7q be retained relative to k 2 v\/ A 2 in order to avoid 
a singularity at k = 0. This means that (z) is localized 
over a larger scale distance than Ao, as indicated earlier. 

The higher cigenfunction 



A k = fcexp (-fc 2 /2) 



(50) 



remains unstable even for Cj. < and corresponds to the 
first unstable mode analyzed in Ref. [5] for Cj. = 0. In 



this case we obtain Tn = 2 



3C£, i.e. 



7o 2 



6vm_ 
7 Ho 



V3 



and 



El 



(3r? T - 2) Q k + 



Ho 



(51) 



(52) 



We observe that the relevant growth rate vanishes for 
Cq < —2/3, that is for relatively flat temperature profiles. 
Moreover, in this case the 7gQ contribution to Eq. (43) 

is not needed to avoid a singularity at k — and Q; oc 
cxp(-P/2). 

VII. Flatter Temperature Profiles. When ijt < 
2/3 and Cq is negative the higher even modes have to be 



considered. For this we define Cg 



-Cq and we note that 



the first relevant solution in this class has the expression 

A k = A° k (1 + a k k 2 ) exp (-P/2) , (53) 

where a necessary condition, in order to find Tq > 0, is 
that a k > 1/2 as can be seen from the asymptotic solution 
of Eq. (41) for k 2 -> 0. 

In particular, by inserting the expression (53) in Eq. 
(41) we obtain 



= C° (2a k - 1) . 



and 



E° k = 7 - C ° = 3 - C ° + a k (T 2 Q + 5C ° - 6) 



Thus 



1 /Tg \ _ 4 
ak - 2 [c% + 1 ) ~ (T 2 + 5C°-6) 



and we obtain the dispersion relation 

r^-6(i-c °)r 2 -c (i4-5C °) = o 



(54) 

(55) 
(56) 

(57) 



3 (1 - C° ) + v/9 (1 - C°y + C ° (14 - 5C °). (58) 



From this we see that for C — > 0, Tq — > 6 and that 
Tq = for Cq = 14/5, that is T decreases monotonically 
until it vanishes as Cq increases from zero to 14/5. When 
C§. < —14/5 a higher eigensolution of Eq. (41) has to 
be considered. We note that, when Cg — > 0, a k — > 3/Cg 
while the corresponding expression for A z (z) remains a 
well behaved function. 

Moreover, we can verify that if Cj. > this mode re- 
mains unstable and that Tq increases as C^p increases. 

VIII. Vertical Transport. — The considered mode 
can produce a significant rate of particle density transport, 
in the vertical direction, that is of contrary sign to that 
of the temperature transport and modify the density and 
temperature profiles in such a way as to lead rjx toward 
2/3. Thus, if tjt > 2/3 a particle inflow toward the equa- 
torial plane is induced while the temperature rises away 
from the equatorial plane. In fact, this is somewhat similar 
to the theoretical explanation [10] for the particle inflow 
observed in magnetically confined toroidal plasmas that is 
associated with the radial electron temperature gradient. 
In the present case we may argue that tjt can be kept 
above 2/3 by a plasma heating process [3], such as one 
due to an anomalous viscosity, that is localized around the 
equatorial plane and attracts particles toward the equato- 
rial plane. A significant cooling process at the edge of the 
plasma disk can reinforce the described type of transport. 

When r/T < 2/3, including the case where t\t — or 
where the surface of the disk can be hotter than the inte- 
rior, a rise of the central temperature is induced while the 
particle transport is outward, from the equatorial plane. 
These arguments arc based on the relevant quasi-linear 
analysis that gives the vertical particle flux 



T pz = ((nv a » ~ --7o 



6 



dz n 2Tdz 



(59) 



where (()) indicates an average over a radial distance AR 
such that 1/fco < AR < Rq. The corresponding tempera- 
ture flux is /(fv z ^\ ~ - ((nv*» T/n. 

IX. Conclusions. — The vertical density and tem- 
perature gradients in a differentially rotating plasma disk 
combined with those of the radial gradient of the rotation 
frequency are the key factors involved in the excitation 
of localized (ballooning) modes that are contained within 
the disk when a seed vertical magnetic field is present. 
The most important effect of these modes is that of in- 
ducing a change of both the density and the temperature 
profiles around the equatorial plane until these reach the 
"adiabatic condition" {dlnT / dz / (dlnn / dz) = 2/3 . Thus 



transport in the vertical direction is shown to have a rele- 
vant role in the dynamics of disks while, previously, only 
the radial fluxes of angular momentum and particles had 
been paid attention to. We note that the class of modes 
considered here benefits from the relatively strong vertical 
gradients of particle density and temperature that can be 
present in thin disks. On the contrary, the well known 
MRI instability that is well suited to the geometry of a 
differentially rotating plasma cylinder is penalized by the 
thin disk geometry as this forces the instability driving 
factory to compete with a relatively strong bending of the 
magnetic filed lines. Finally, the class of modes that we 
have analyzed can lead to the formation of a non-linear 
magnetic "crystal" structure, within the disk, made of a 
sequence of counter streaming toroidal current channels 
as described in Rcf. [11] and may start the evolution of 
a classical disk configuration toward a sequence of plasma 
rings [4]. 
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